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Abstract 

Exact solution of the Schrodinger equation with deformed ring shaped potential 
is obtained in the parabolic and spherical coordinates. The Nikiforov-Uvarov method 
is used in the solution. Eigenfunctions and corresponding energy eigenvalues are cal- 
culated analytically. The agreement of our results is good. 
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1 Introduction 



The exact solutions of the Schrodinger equation with some well-known central and noncen- 
tral potentials is an important subject in quantum mechanical problems. Such solutions are 
helpful for checking and improving models and numerical methods besides of understanding 
about the characteristics of a quantum system. The exact solutions of this equation are 
possible only for certain potentials such as Coulomb, Morse, Poschl- Teller and harmonic 
oscillator etc. [1]. The other exactly solvable one is the deformed ring-shaped potential 
introduced by Hartmann [2]. This potential involves an attractive Coulomb potential with 
a repulsive inverse square potential one. In spherical coordinates it can be defined as 



V{r,6) 



2 a 
- q5- 



5<7 2 a e , (1) 
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where ao and eo denote the Bohr radius and the ground state energy of the hydrogen atom 
respectively. 5 and a are positive real parameters as well. Their range varies from 1 upto 
10. This potential can be used in quantum chemistry and nuclear physics to describe ring- 
shaped molecules like benzene and interactions between deformed pair of nuclei. We point 
out that the potential takes the form of the Coulomb potential in the limiting case 5a 2 = Z 
and q = for hydrogen-like atoms. The energy eigenvalues of the potential has been 
calculated before by using some useful methods. For example, these are a non-bijective 
canonical transformation, namely, Kustaanheimo-Stiefel (KS) transformation, dynamical 
group method, path integral and SUSYQM method etc. [2-17]. Moreover, this potential 
can be defined as a Coulomb plus Aharonov-Bohm potential by defining the parameters as 
-eef = 2a e 5a 2 , -A/2/i = qe a 2 5 2 a 2 and B = 0. [18, 19, 20]. 

In the present work, we have obtained an exact solution of the Schrodinger equation with 
the q-deformed ring shaped potential by using the Nikiforov-Uvarov (NU) method in both 
parabolic and spherical coordinates. The method is based on the solution of differential 
equation transformed into the hypergeometric type [21, 22]. 

The paper is organized as follows: In Section II we introduce the Nikiforov-Uvarov 
method. In Section III we apply the method to solve the Schrodinger equation in both 
parabolic and spherical coordinates respectively. In Section IV we present numerical results 
for Z = 1 and q = with the conclusion. 



2 The Nikiforov-Uvarov Method 



The Nikiforov-Uvarov method first reduces the second order differential equations (ODEs) 
to the hypergeometric type with an appropriate coordinate transformation x = x(s) as 

ds 2 + a(s) ds ^ a 2 (s) {) { ) 

where a(s) and a(s) are polynomials with at most second degree, and f(s) is a polynomial 
with at most first degree [21, 22]. If we take the following factorization 

= <f>(s) y(s), (3) 

the Eq. (2) takes the form [22] 

^)^ + t( s) ^W)=0, (4) 



where 



and 



( x( S )= 7 r( S )^(ln0( S )), (5) 



T (s) = f(s) + 2?r(s). (6) 



Also, A is given 

A„W + '" ( "- 2 ^"' =0, „ =0ili2 ,... {7) 

The energy eigenvalues can be calculated from the above equation. We first have to 
determine n(s) and A by defining 

fc = A-7r'(s). (8) 
Solving the quadratic equation for n(s) with the Eq. (8), we get 



a - t \ • l a' — r x 2 



Here, 7r(s) is a polynomial with the parameter s and prime factors denote the differentials 
at first degree. The determination of k is the essential point in the calculation of n(s). It is 



obtained by setting the discriminant of the square root to zero [22]. Therefore, we obtain 
a general quadratic equation for k. 

The determination of the wave function is now in order. We consider the Eq. (5) and 
the Rodrigues relation 

*M = ^£ [•"«*>]• < 10 > 

where C n is normalization constant and the weight function p(s) satisfy the following rela- 
tion 

±[a(s)p(s)}=T(s)p(s). (11) 

The Eq. (10) refers to the classical orthogonal polynomials that have many important 
properties especially orthogonality relation can be defined as 

rb 

/ Vn(s) y m (s)p(s) ds = 0, m^n. (12) 

J a 

3 Calculations 

The Schrodinger equation in spherical coordinates becomes 

2th 

V 2 ^ + [E - V{r,6)} = 0. (13) 

h 

We are first going to study for solution of the problem using parabolic coordinates. 
3.1 Parabolic coordinates 

One can write the second type parabolic coordinates as [3, 4, 23] 

x = costp, y = £rj sirup and z = ]- (if - £ 2 ) , (14) 



and 



£ 7] — r sin9 and r = — (if + £ 2 ) . (15) 
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If we write trial wave function in the following form [4] 



^^, v , V ) = -^=u(0v( V )e^, (16) 



one can get two-coupled differential equations 



d 2 u (T 2 -i) ,M U2 ,2m. 



and 



cr 2 *; (T 2 -i)^ { ,2m 



u + (-£2-)»7 « - = 0, (18) 



drf if h h 



where T = \Jm' 2 + qS 2 a 2 and — 2<r 2 <5 ao eo- We will first solve the Eq. (17) and 

then easily get the other one. 

By using the transformation s = £ 2 , the Eq. (17) is therefore transformed into the equation 
of hypergeometric type. Hence, we have 

u"(s) + —u'(s) + At \-e 2 s 2 - a 2 s - /3 2 ] u(s) = 0, (19) 
2s 4s 2 L J 

2 2mE 2 2m A a 2 1^2 ^ 
where £ = — 3—, a 1 = — rp^i an d P — (I )• 

Comparing the Eq. (19) with the Eq. (2), we get 

a(s) = 2s, f(s) = 1 and a(s) = (-e 2 s 2 - a\s - (3 2 ) . (20) 
Substituting these into the Eq. (9), we write 



tt(s) = -± -^Ae 2 s 2 + (8k + 4a 2 ) s + A(3 2 . (21) 



The constant k can be determined from the condition that the discriminant of the square 
root must be zero, so that 

k h2 = ~a 2 ±ep. (22) 
Hence the final result for Eq. (21) can be written as 



, ( (es — 0), for k = —\a\ — ef3 
7r(*) = -± " (23) 
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(es + /3), /or k = -\a{ + el3. 



A proper value for tt(s) is chosen, so that the function 



t(s) = 2(1 + (5) - 2es, 
must have a negative derivative [22]. From the Eq. (7) we can obtain 



(24) 



An = ~0l\-£- @E 

= 2ne. (25) 
Following the same procedure again one gets for Eq. (18) as 



A„/ — 



-—a n 
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e-f3e 



(26) 



By combining each side of the Eqs. (25) and (26) we obtain energy eigenvalues 



E„ 



S 2 a 4 



e . 



n,ri = 0,1,2 



(27) 



_(n + n' + l + py 

This solution is identical for (3 ~ T 2 with the ones obtained before [3, 4, 6, 7, 8, 9, 13]. 

Now, we are going to determine the wave function. Considering the Eq. (3) and using 
the Eq. (5) we get 



^(s) = s u/4 e-i s , (28) 

where v — 1 + 2(3. 

From the equations Eqs. (11) and (10), we obtain 

*" (s > = ^ (29 > 

with p(s) = s^ -1 )/ 2 e~ £S . The Eq. (33) stands for the associated Laguerre polynomials. 
That is 



Vn{s) = L?(s), 

Hence we have found the wave function that belongs to the Eq. (17) as 



(30) 



with s = £ 2 . Similarly, we can also write the wave function for Eq. (18) 



(31) 



v n (v) = C n . a"' 4 e-f (32) 
with s — rj 2 . Therefore, the total wave function takes 

*n,n>Mt,V,<P) = ~k= C n,n> e~" L^s) L^s) (33) 



where the normalization constant C n>n i can be found from the Eq. (12) as 



4 « ,,,/ 0,1,2... (34) 



\ (n + (3)\ (n' + (3)V 

One can easily see that in the case of rsinO = £r), the problem reduces to harmonic 
oscillator plus inverse square potential case. The latter, we have studied that this problem 
also reduces to that of molecular Kratzer potential like (Coulomb plus inverse square). 

3.2 Spherical coordinates 

Considering the Eq. (13) we write the total wave function as 

f(r,^) = ^e(0)%), (35) 
with the well-known azimuthal angle solution 

Sfa) = -L m = 0,±l,±2,... (36) 

Thus we get as 



1 d ( JQ\ ( (m 2 + b')\ „ 

•sine— ) + [k- K - — i 8 = 0, (37) 



and 
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d 2 U 2 / . a' k \ tt 

+ -1-E' o)U = 0, (38) 



dr 2 7 \ r r 2 



2mE . 2ma . 2mb , 

where h = — ^- , a = — ^- and b = — k and m are also seperation constants. 
n n n 



Using the NU-method, we are going to solve them. By defining m! = \/m 2 + b in Eq. 
(37) and taking x = cos9, it will have a form of hypergeometric type 



d 2 B 2x dQ 



+ 



dx 2 (1 — x 2 ) dx (1 — x 2 ) 2 
Comparing it with the Eq. (2), we get 



k{\ -x 2 ) -m' 2 G(x) = 



a(x) = x, f(x) = —2x and a(x) = k(1 — x 2 ) — m' 2 . 
Substituting these into the Eq. (9), we get 



(39) 



(40) 



ir(x) = ±yf-(k + «)(1 -x 2 ) +m' 2 . (41) 

The constant k is determined from the condition that the discriminant of the square root 
must be zero. Thus, we find 



± m', for k = k 

ir(x) = < 

± m' x, for k = k — m' 2 . 
A proper value for n(x) can be chosen, so that the function 

T (x) = -2(m' + 1) x, 
has a negative derivative. From the Eq. (7) we can obtain 



(42) 



(43) 



\ n = k — m'{m' + 1) 

= 2n (m' + l) + n (n - 1). (44) 

Solving for k, we have 

K = K n = £'(£' + 1), (45) 

where £' — n + m f . 

Now we are going to determine the wave function. From the Eqs. (11) and (10), we can 
write 



y n (x) = 



(1-x 



l\n+m' 



with p(x) 



p(x) dx n 

[1 — x 2 ) m ' . Therefore the Eq. (46) stands for Jacobi polynomial as 



(46) 



where n = I' — m' . The wave function becomes 



(47) 



Q(x) = Ql>, m > 

= c v , m> {i-xY' 2 Pt-S\^ 

with x = cos9 (x G [—1, 1]). Using the Eq. (12), we get 



Ci'm' — 



I2£' + 1 



2 m '(f + l)V 2 
Let us now consider the Eq. (38) 



(£' - m')\(£' + m')\ 



(48) 



(49) 



u"(r) + 



—E'r 2 — a'r — k 



u(r) = 0, 



Comparing this equation with the Eq. (2), we obtain 



a(r) = r, r(r) = and a{r) = —E'r 2 — a'r — k. 
If we insert these into the Eq. (9), one gets 



?r(r) = \ ± ^AE'r 2 + Mk + a')r + (1 + 4k). 



(50) 



(51) 



(52) 



We can determine the constant k by using the condition that discriminant of the square 
root is zero, that is 



fci, 2 = -a'±yjE' (1 + 4k). 
Hence the final form of the Eq. (52) for each value of k becomes 



(53) 



1 l [2y/Wr + VI + 4k], /or fc = -a' + y/E'(l + 4k) 
[2-/BV- V1 + 4k], /or fc = -a' - ^E'(\ + Ak). 



(54) 



A proper value for n(r) is taken, so that the function 



r(r) = (1 + Vl + 4k) - 2V£V, 
has a negative derivative. From the Eq. (7), we can write 



(55) 



A n = -a' - /E'(1+4k)- v 7 ^ 7 
= 2n v 7 ^ 7 . 



(56) 



Therefore, it gives us the energy eigenvalues of the radial equation with the deformed ring- 
shaped potential 



E 



S 2 a 4 



n r + £' + 1 



, n r = 0,1,2,... 



(57) 



where n r denotes the radial quantum number which belongs to the Eq. (38). 

To determine the wave function, we consider the Eqs. (3) and (5) for obtaining 



0(r) = e 



E> r „(i/-l)/2 



(58) 



where v — 1 + 2 a/1 + 4k. Thus from the Eqs. (11) and (10) we have 



Vn{r) 



a n{r) p(r) 



(59) 



p(r) dr n 

with p(r) = z~^® r r^~ v ^ 2 . The Eq. (59) stands for associated Laguerre polynomials, that 
is 



y n (r) = L k n (r), 

where k — {y — l)/2. The radial part wave function is written as 

U n (r) = C n e-^ r' k L*(r). 
By using the orthogonality condition, we can determine the coefficient as 



a, 



TV. 



yn ' k y 2(n + k)(n + k)V 
with k — £'(£' + 1). Hence, the total wave function takes the form 



(60) 



(61) 



(62) 



9(r,6,<p) = 4= Cn,k C t ,, m ,\e-^ r k (Sinfi)' 
V Ztt l 



x p( m '' m '\cos9) L*(r) e im ^] . (63) 

4 Conclusion and Remarks 

We have obtained the exact eigenfunctions and corresponding energy eigenvalues of the 
Schrodinger equation with the deformed ring-shaped potential in both second type parabolic 
and also spherical coordinates by using the Nikiforov-Uvarov method. At first our problem 
reduces to the harmonic oscillator plus inverse square potential, it also reduces to the 
problem that molecular Kratzer (Coulomb plus inverse square) one in second case. Results 
obtained in two different coordinate systems are identical by following the conditions (3 ~ T 2 
in Eq. (27) and £' — n+m' in Eq. (57). Some numerical values of energy for a hydrogen-like 
atom due to the attractive Coulomb potential are presented in tabular form. The total wave 
functions, in both coordinates, are physical. They behaves asymptotically. The agreement 
of our analytic and numerical results is good. 
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Table I: Energy levels for a hydrogen-like atom with 5a 2 = Z and q = 0. 
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